Centaurus A as the Source of ultra-high energy cosmic rays? 



Claudia Isola", Martin Lemoine'', Giinter Sigl^ 

" Centre de Physique Theorique, 
Ecole Polytechnique, 91128 Palaiseau Cedex, France 
' Institut d'Astrophysique de Paris, C.N.R.S., 98 his boulevard Arago, F-75014 Paris, France 

(April 13, 2001) 



O 
O 

< 



> 

oo 
(N 

O 



Oh. 

6 ■ 



We present numerical simulations for energy spectra and angular distributions of nucleons above 
lO'^^ eV injected by the radio-galaxy Centaurus A at a distance 3.4 Mpc and propagating in extra- 
galactic magnetic fields in the sub-micro Gauss range. We show that field strengths B ~ O.SfiG, as 
proposed by Farrar & Piran, cannot provide sufficient angular deflection to explain the observational 
data. A magnetic field of intensity B ~ IfiG could reproduce the observed large-scale isotropy and 
could marginally explain the observed energy spectrum. However, it would not readily account for 
the E = 320±93EeV Fly's Eye event that was detected at an angle 136° away from Cen-A. Such a 
strong magnetic field also saturates observational upper limits from Faraday rotation observations 
and X-ray bremsstrahlung emission from the ambient gas (assuming equipartition of energy). This 
scenario may already be tested by improving magnetic field limits with existing instruments. We 
also show that high energy cosmic ray experiments now under construction will be able to detect the 
level of anisotropy predicted by this scenario. We conclude that for magnetic fields B ~ 0.1 — 0.5 fiG, 
considered as more reasonable for the local Supercluster environment, in all likelihood at least a few 
sources within ~ 10 Mpc from the Earth should contribute to the observed ultra high energy cosmic 
ray flux. 



I. INTRODUCTION 

In acceleration scenarios ultra high energy cosmic rays 
(UHECRs) with energies above 10^® eV are assumed to 
be protons accelerated in powerful astrophysical sources. 
During their propagation, for energies above ^ 50 EeV 
{lEeV = lO^^eV) they lose energy by pion production 
and pair production (protons only) on the microwave 
background. For sources further away than a few dozen 
Mpc this would predict a break in the cosmic ray flux 
known as Greisen-Zatsepin-Kuzmin (GZK) cutoff [1], 
around 50 EeV. This break has not been observed by 
experiments such as Fly's Eye [2], Haverah Park [3], 
Yakutsk [4] , Hires [5] and AGAS A [6] , which instead show 
an extension beyond the expected GZK cutoff and events 
above 100 EeV. This situation has in recent years trig- 
gered many theoretical explanations ranging from con- 
ventional acceleration to new physics as well as the con- 
struction of large new detectors [7] . 

In bottom-up scenarios of UHECR origin, in which 
protons are accelerated in powerful astrophysical objects 
such as hot spots of radio galaxies and active galactic nu- 
clei [8] , one would expect to see the source in the direction 
of arrival of UHECRs. The lack of observed counterparts 
to the highest energy events [9,10] implies the existence 
of large scale intervening magnetic fields with intensity 
B ~ 0.1 — 1 iiG [10], which would provide sufficient an- 
gular deflection, or bursting sources and a magnetic fleld 
of intensity B > 10~"'^^G which would impart sufficient 
time delay to UHE protons to explain their lack of cor- 
relation in time of arrival with optical or high energy 



photons [11]. It has been realized recently that mag- 
netic fields as strong as ~ 1/x G in sheets and filaments 
of large scale structure, such as our Local Supercluster, 
are compatible with existing upper limits on Faraday ro- 
tation [12-14]. 

Such strong magnetic fields would have profound con- 
sequences for the propagation of charged ultra-high en- 
ergy cosmic rays [15-17]. In particular in the presence 
of a magnetic fleld of strength B ~ 0.1 /iG, and for 
a source at distance ~ 10 Mpc, UHECRs with energy 
i? < 100 EeV would diffuse, while higher energy cosmic 
rays would propagate rectilinearly. The resulting modi- 
fication of the energy spectrum would explain naturally 
the observed spectrum for a unique power-law injection 
of index ~ 2.2. In Ref. [17], we further showed that 
a large number of such sources in the Local Superclus- 
ter, at distance scale ~ 10 Mpc and an ambient magnetic 
field i? ~ 0.1 /iG would explain the large scale isotropy 
of arrival directions observed by AGASA [6]. It would 
also explain the observed small-scale angular clustering 
(five doublets and one triplet within 2.5° out of 57 events 
above 40 EeV) by magnetic lensing effects through the 
large scale turbulent magnetic field. 

At first sight this suggests the possibility of having only 
one object in the Sky as the source of all UHECRs with 
E > 5EeV, including the highest energy events. Two ver- 
sions of such single source scenarios have recently been 
put forward in the literature: one with an extreme ver- 
sion of a coherent Galactic magnetic wind structure in 
which all observed UHECRs supposedly can be traced 
back to M87 in the Virgo cluster [18], and a second 
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one with Ccntaurus A at 3.4 Mpc distance with an all- 
pervading magnetic field of intensity B ~ 0.3 /iG [19]. In 
the present paper we examine critically the latter of these 
scenarios using detailed numerical simulations for the en- 
ergy spectrum and the angular distribution of UHECRs 
propagating in magnetic fields of r.m.s. strength between 
0.3 and 1 iJ,G. 



II. NUMERICAL SIMULATIONS 

Trajectories are calculated from the Lorentz equation 
in a given magnetic field and pion production is treated 
as stochastic energy loss while pair production is included 
into the equations of motion as a continuous energy loss 
term. This numerical tool has been used and discussed 
in earlier publications [16,17]. 



nms that it has no effect on the results discussed in the 
following. The magnetic field modes are dialed on a grid 
in momentum space according to this spectrum with ran- 
dom phases and then Fourier transformed onto the cor- 
responding grid in location space. The r.m.s. strength B 
is given by = dkk"^ (^^(A;)). 

Typically, 5000 trajectories are computed for each 
magnetic field realization obtained in this way for 10- 
20 realizations in total. Each trajectory is followed for a 
maximal time of 10 Gyr and as long as the distance from 
the observer is smaller than double the source distance. 
We have also checked that the results do not significantly 
depend on these cut-offs. Furthermore, the distance limit 
is reasonable physically as it mimics a magnetic field con- 
centrated in the large scale structure, with much smaller 
values in the voids, as generally expected. 
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FIG. 1. The angular image in terrestrial coordinates, aver- 
aged over all 20 magnetic field realizations of 5000 trajecto- 
ries each, for events above 40 EeV (upper panel) and above 
100 EoV (lower panel), as seen by a detector covering all 
Earth, for the case suggested in Ref. [19] corresponding to 
B = 0.3 /jG, and the source Cen-A located 3.4 Mpc away. 
The grey scale represents the integral flux per solid angle. 
The solid line marks the Supergalactic piano. The pixel size 
is 1°; the image has been convolved to an angular resolution 
of 2.4° corresponding to AGASA. 

We assume a homogeneous random turbulent magnetic 
field with power spectrum {B{kY) <x. fc"^ for 27r/L < k < 
27r//c and (^^(A:)) = otherwise. We use ns = —11/3, 
corresponding to Kolmogorov turbulence, in which case 
L, the largest eddy size, characterizes the coherence 
length of the magnetic field; we use L ~ 1 Mpc, which 
corresponds to a few turn-arounds in a Hubble time. 
Physically one expects Ic ^ L, but numerical resolution 
limits us to Ic > 0.008L. We generally use Ic ^ 0.03 Mpc, 
but we checked by increasing the resolution for several 
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FIG. 2. The distribution of arrival declination on Earth, 
averaged over many realizations, for E >40 EeV (dotted line) 
and E > 100 EeV (solid line), for the scenario corresponding 
to Figs. 1. The dash-dotted line represents an isotropic dis- 
tribution. The pixel size is 1° and the image has again been 
convolved with an angular resolution of 2.4°. 



In Fig. 1, we show the angular distribution of UHE 

events as seen on Earth in equatorial coordinates, for 
two ranges of energies E > 40 EeV, and E > 100 EeV, 
and for B = 0.3 /iG, with Cen-A as the source. These 
images are averaged over different spatial realizations of 
the magnetic field. The distributions for specific realisa- 
tions arc more anisotropic than the average due to cosmic 
variance. The angular distribution predicted by this one 
source model is thus not consistent with the isotropic 
distribution deduced by experimental data [2,6]. This is 
even more clearly demonstrated by Fig. 2, which gives the 
distribution in declination of arrival coordinates of UHE 
events. The source Cen-A is located at RA = 201.3°, 
6 = —43.0° in equatorial coordinates, and corresponds 
to the peak of flux in Figs. 1,2. 
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FIG. 3. The realization averaged energy spectra corre- 
sponding to Figs. 1,2. The solid line represents the spectrum 
that would have been detected by AGASA, and has been ob- 
tained from Eq. (1). The dashed line indicates the spectrum 
uniformly averaged over the whole sky. The dotted line is the 
spectrum predicted by an AGASA type experiment in the 
Southern hemisphere. The one sigma error bars indicate the 
AGASA data. The various spectra have been normalized to 
optimally fit the AGASA flux. Significant uncertainties due 
to cosmic variance and parameters such as the largest eddy 
size L only occur for energies between ~ 70EeV and a few 
hundred EeV where they are still smaller than a factor ~ 2.5. 



production of higher energy particles, also contributes to 
the change of slope at low energies. 

The scenario just discussed, which is already ruled out 
by the energy spectrum and large-scale isotropy recorded 
by Northern hemisphere detectors, has been proposed by 
Farrar & Piran [19] on analytical grounds to explain all 
observational data. The discrepancy, as will be discussed 
in greater detail in Section 3, results from the fact that 
they argued that diffusion held up to the highest ener- 
gies, whereas in fact the diffusion approximation breaks 
for E > 100 EeV, implying much larger anisotropics at 
these energies. The impact of angular anisotropy on the 
energy spectrum for a source located in the blind area of 
Northern hemisphere detectors had also been overlooked 
in Ref. [19]. 

For the scenario with B = 0.3/z G we have also estab- 
lished a rough estimate of the minimal number of sources 
necessary to explain existing observations in the follow- 
ing way: We overlaycd distributions from single sources 
in random directions (for simplicity chosen to be at equal 
distances as Cen-A) and computed the average and fluc- 
tuations of the declination distributions of UHECR ar- 
rival directions. The minimal number of sources is de- 
termined by requiring that the distribution is consistent 
with isotropy within the fluctuations and turns out to be 
5-10. 



Cen-A does not lie in the field of view of experiments 
that have provided data so far such as the Fly's Eye and 
AGASA. This and the fact that the angular deflection of 
particles with E > 100 EeV is relatively small has an im- 
portant consequence for the energy spectrum predicted 
by this model: The Northern hemisphere experiments 
should never have detected the highest energy events, 
for which the angular deflection is too weak to bring the 
particle in the field of view. This is made clear in Fig. 3, 
where we compare the spectrum predicted for AGASA 
with the actually observed spectrum, assuming an injec- 
tion spectrum cx E^"^ up to 10^^ eV. The prediction is 
obtained by folding the simulated distribution in energy 
E and arrival direction 17, D{E, 17), with the normalized 
AGASA exposure function AGASA((5) (which to a good 
approximation only depends on declination 5): 

j{E) = J dnD{E,6)AGASA{S). (1) 

For reference, we also show in Fig. 3 the spectrum that 
would be observed by an idealized detector covering the 
whole sky uniformly, and by a mirror AGASA experi- 
ment located in the Southern hemisphere. The solid an- 
gle integrated spectrum J dflD{E,S) observed by a uni- 
form detector is still different from the injection spectrum 
(oc E~^) at low energies, while the two match at high en- 
ergies. This results from an increased local residence time 
due to diffusion at low energies, and rectilinear propaga- 
tion (hence unaffected energy spectrum) at high ener- 
gies [16]. The pile- up around E ~ 40 EeV due to pion 
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FIG. 4. The distribution of arrival declination on Earth, 
averaged over many realizations, for E > 100 EeV (dotted 
line) and E > 300 EeV (solid line) . All other parameters are 
as in Fig. 2, expect B = 1 /uG. 

Wc now investigate whether stronger magnetic fields, 
by providing larger angular deflection, might provide a 
better match to the observational data. In particular, we 
focus on the case where B = IfiG, and Cen-A is again 
the unique source of UHECRs. The resulting arrival dis- 
tribution in declination is shown for several ranges of 
energies in Fig. 4, and the resulting energy spectra, cal- 
culated as before, is shown in Fig. 5. Increasing the 
magnetic fleld strength increases the maximal energy 
at which diffusion takes place, hence it decreases the 
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anisotropy at each energy up to that maximal energy, 
and thus reduces the difFerenccs between the spectra seen 
by different detectors (AGASA, uniform, and Southern 
hemisphere analogue of AGASA). 
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FIG. 5. Normalization averaged energy spectra with the 
same conventions as in Fig. 3, but for i? = 1 fiG, all other 
parameters being equal. 

The large-scale anisotropy in this case is much smaller 
and could not have been detected by Northern hemi- 
sphere experiments such as the Fly's Eye and AGASA. 
The predicted energy spectrum for AGASA does provide 
a very good fit to the observed spectrum, see Fig. 5, but 
the difference is not sufficient to rule out this scenario on 
this ground. However, we note that it is unlikely that the 
highest energy Fly's Eye event would have been detected 
in this model. This event of energy £■ = 320 ± 93 EeV 
makes an angle of 136° with the line of sight to Cen- 
A [10]. By folding the energy probability distribution for 
the Fly's Eye event with the simulated distribution of de- 
flection angles and energies one finds that 95% of events 
with a recorded energy corresponding to the Fly's Eye 
event are deflected less than ~ 130°. 

Thus the present observational evidence is not suffi- 
cient to rule out with a high degree of confidence the pos- 
sibility that Cen-A is the source of all ultra-high energy 
cosmic rays, if the intervening magnetic field B > 1 /uG. 
Notably, one cannot rule out the possibility of having a 
magnetic configuration different from Kolmogorov tur- 
bulence, in which case the modification of the energy 
spectrum would be different than what is shown in the 
present simulations. For instance, it has been proposed 
that the magnetic field is not all pervading, but strongly 
enhanced in regions close to radio-galaxies that were ac- 
tive in the past [10,20], in which case the configuration 
seen by UHECRs would be a collection of scattering cen- 
ters rather than Kolmogorov turbulence. 

However, if future or ongoing experiments in the 
Northern hemisphere, e.g. the High Resolution Fly's 
Eye [5] and AGASA [6] keep recording cosmic rays above 
~ 200 EeV with large deflection angles from the line 



of sight to Cen-A, even the scenario with B ~ 1 fiG 
would be unequivocally ruled out. Similarly, if no sig- 
nificant anisotropy is seen between these experiments 
and the Pierre Auger project [21] at the highest en- 
ergies, the model would be discarded. As an exam- 
ple, we calculate from Eq. (1) the fractional difference 
51 = (In — Is) /{In + Is) of integral fluxes In and Is 
that would be seen by detectors in the Northern and 
Southern hemispheres above a given energy (for simplic- 
ity, we use exposure functions for AGASA and an anal- 
ogous one for the Southern hemisphere, as before). We 
find SI ~ -0.19 for E > 100 EeV, and SI ~ -0.78 for 
E ^ 300 EeV. These numbers can also very roughly been 
estimated from Fig. 4. Anisotropics of this size should 
be easily detectable by a full sky observatory such as the 
Pierre Auger project [21] which is currently under con- 
struction. 

We note in this context that wc have assumed in all 
our simulations that the source emits only protons. If the 
source is sufficiently compact, protons could convert into 
neutrons within the source. As pointed out in Ref. [22], 
for sources as close as Cen-A, neutrons at the highest en- 
ergies could survive decay and produce a spike in the di- 
rection of the source. This can only increase anisotropy. 
Preliminary simulations performed with our code indi- 
cate that the total flux in the ~ 2° pixel centered on the 
source is significantly (i.e. by a factor ~ 2) increased 
only above ~ 300 EeV. 

We stress that S ~ 1 /xG corresponds to the upper 
limit inferred on the strength of the magnetic field in 
the Local Supercluster from Faraday rotation observa- 
tions of distant quasars [12-14]. Strictly speaking, the 
rotation measure is a function of B^/L, where L denotes 
as before the coherence scale of the field, provided that 
L <^ R, where R represents the size of the medium per- 
vaded by the magnetic field (for the Local Supercluster 
i? c± 10 — 50Mpc). In principle, the coherence length L 
could be smaller than 1 Mpc, used in the previous Sec- 
tion, and B correspondingly larger. However, as will be 
discussed in the next Section, the diffusion coefficient 
D{E) scales inversely proportional with L, so that by 
decreasing L, one would increase D{E) correspondingly, 
and the diffusion approximation would break down at 
a smaller energy, thereby increasing the anisotropy at 
higher energies. We thus find that one cannot decrease 
L and increase B to improve the fit to the data. More- 
over, notwithstanding this fact, if diffusion were to be 
more efficient at the highest energies observed, the effec- 
tive distance traveled would be consequently increased, 
and the GZK cut-off would be more pronounced, which 
would further aggravate the disagreement with the ob- 
served spectrum. 

We thus conclude that B ~ 1 /xG seems to be the only 
value of B that would be marginally consistent with Cen- 
A as the source of observed UHECRs. Wc further note 
that if equipartition of energy holds, the magnetic field 
intensity is tied to the thermal energy density of the am- 
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bient gas: 

B = 0.5nGT^^^KYQ^h7Q, (2) 

where T7 = T/10''K is the temperature of the Local Su- 

pcrclustcr in units of 10^ K, and kiq = k/10, with k 
the collapse factor (i.e., the local overdensity of baryons 
and electrons), and /170 the Hubble constant in units of 
70 km/s/Mpc. The gas cools by brcmsstrahlung emission 
in the kcV range, which in principle can be observed. 
A marginal detection of X-ray emission correlated with 
the plane of the Local Supercluster has actually been 
reported and corresponds to a collapse factor kiq ~ 1, 
with a weak dependence on the assumed temperature 
T ~ 10^ K [23]; the signal is however weak and these 
parameters could be in error. Numerical simulations of 
large-scale structure formation indicate that kio — 1 and 
T ~ 10 ''K are probably upper limits for sheets such as 
the Local Supercluster, and seem to better describe the 
filamentary structures [13]. Deep searches for soft X-ray 
emission correlated with the plane of the Local Super- 
cluster using the XMM-Newton or Chandra observatories 
could improve these limits, and appear mandatory. 

III. COMPARISON WITH ANALYTICAL 
ESTIMATES 

Let us now compare these results with analytical es- 
timates in an approach similar to [19]. There is often 
confusion in the literature about different regimes of dif- 
fusion and corresponding expressions for the diffusion co- 
efficient. It is dangerous to take analj^tical expressions 
too literally as there exists no analytical derivation of 
diffusion coefficients in the limit in which the turbulent 
component of the magnetic fields becomes comparable 
or stronger than a (putative) uniform component, the 
so-called strong turbulence regime. In this respect, one 
should note that the formula for the diffusion cocfhcicnt 
given in Ref. [15], now used without caution in the com- 
munity, does not correspond to an analytical derivation. 
It is a phcnomcnological formula, that is furthermore 
shown to be in error in Ref. [24], where analytical ap- 
proximations and accurate measurements of diffusion co- 
efficients obtained through Monte-Carlo simulations are 
presented. Notably, this latter study shows that for Kol- 
mogorov turbulence the diffusion coefficient can be ap- 
proximated as: 

D{E) ~ 0.02E2qBzI MpcVMyr E20 < E^, 

~ QmEl^BzlL-'^ MpcVMyr E20 > E^. (3) 

In this expression, £20 is the UHECR energy in imits of 
100 EeV, B-Q is the magnetic field strength in units of 
HG, L is in units of Mpc and Ec = 1A5B-qL (in units of 
100 EeV) corresponds to the condition tl = L/2tt, where 
tl = 0.11i52o-BZ6 Mpc is the Larmor radius. Note the 
difference of the above result with the formula given in 



Ref. [15], for which D{E) oc E^/^ for £^20 < E^, and 
D{E) oc E for E20 > E^. The dependence of D{E) for 
E < Ec in Eq. (3) above agrees with the phenomenolog- 
ical Bohm diffusion coefficient Db — r^c within a factor 
two. 

The diffusive regime as well as the transition to nearly 
rectilinear propagation can also be seen in the depen- 
dence of time delay r (defined as the difference between 
the total propagation time and the straight flight distance 
d/c) on energy E shown in Fig. 6 for B_q = 1 and Cen-A 
as the source: In the diffusive regime the average time 
delay t{E) ~ (P /AD{E), where d is the source distance, 
whereas for £^20 ^ 3, in the regime of almost rectilinear 
propagation, t{E) ~ 1.7 x 10Vr^2"o'4cnA^S-6 
where dcenA = (i/3.4Mpc. The values obtained for 
D{E) ~ d'^/4:T{E) from Fig. 6 in the diffusive regime 
are consistent with Eq. (3) within the width of the dis- 
tribution. 




ergy E, for B = 1 /iG, and Cen-A as the source. 

Note that according to Eq. (3) a largest eddy size L 
considerably smaller than 1 Mpc (which is not excluded) 
would lead to less diffusion, tending to make anisotropics 
even larger. 

Using these results, one can understand the situation 
encountered in the previous section. In the diffusive 
regime one finds that the ratio of the diffusive propa- 
gation time (which is equal to the time delay in this 
approximation) to the source distance reads: T{E)/d ~ 
16i?^^i?^gLdcenA, for E2Q > Ec- Diffusion ceases to be 
a good approximation when T{E)/d is no longer » 1. 
For instance, T{E)/d > 2 requires £20 < 3i3_6£^/^. 
This corresponds to a maximal energy for diffusion of 
~ 100 EeV when B ~ 0.3 /xG, and ~ 300 EeV when 
B a 1 jiG. These numbers are in agreement with the 
results of the previous section. Note that the ratio of the 
diffusive propagation time to the source distance also cor- 
responds to the ratio of the source distance to the mean 
free path for scattering on the magnetic inhomogeneities 
~ D/c within a factor four. 



5 



The difFcrcncc bctwen our results and Ref. [19] can be 
explained as a misuse of the diffusion coefficient given 
in Ref. [15] by the authors of Ref. [19], and to some 
extent, by the fact that that diffusion coefficient itself 
has the wrong scaling with E, B, and L. More pre- 
cisely in Ref. [19] the authors use the diffusion coefficient 

D{E) ~ OmUpc^/MyrEll^^BZl^^L^/^ [15] (note that 
the present prefactor is the correct one; their Eq. (1) has 
a numerical error) on the whole energy range, whereas 
according to Ref. [15], this scaling is only valid when 
tl < L/2tt, or as above E20 < E^ = 1.45S_6-£/. This 
makes an important difference, because, had the au- 
thors of Ref. [19] used the other limiting regime they 
quote, namely D ~ 0.1E2oBzlMpc^ /Myi, valid for 
£^20 ^ lA5B-eL, they would have found that the ratio 
of the diffusive distance traveled to source distance reads 
d/AD ~ 2.6£'^^B_6rfcenA, which shows that S > 1 /xG is 
necessary to achieve diffusion up to the highest energies 
E20 — 3. Instead, they used the former approximation, 

giving d/AD ~ b.SE^o^^Bl^^L-y^dcenA, which due to 
the (incorrect) weak scaling would let believe that diffu- 
sion can be achieved easily with B ~ 0.3 /xG on the whole 
range of energies. 

One should mention that the transition between dif- 
fusive and rectilinear regimes of propagation is not sud- 
den; it stretches over half an order to an order of magni- 
tude. As one approaches this transition, the anisotropy 
increases steeply to match the small angle deflection in 
the rectihncar regime 0e ^ U0°E2q B^eL^^'^d^^^j^ [11]. 
Therefore it is incorrect to assume that the diffusive esti- 
mate for the anisotropy remains valid up to the transition 
energy. 

Finally, Ref. [19] (unlike Ref. [22]) also base their cal- 
culations for the flux of particles detected at Earth on 
a time-dependent solution of the diffusion equation, but 
one should rather use a stationary solution corresponding 
to a continuously emitting source. Indeed if Cen-A were 
a bursting source, or more generally a source emitting 
only once on a timescale fem Tminj where Tmin is the 
smallest time delay imparted to all UHECRs, or equiv- 
alently, the time delay at the highest energies observed, 
an experiment like AGASA would record events only in 
a limited range of energies [25] . In other words the dis- 
tance do of the diffusion front from the source at time 
t is given by dn{E) = 6D{E)t, and depends on energy. 
At the present time at which AGASA is operating, the 
front of low energy particles would not have yet reached 
the Earth, while that of higher energy particles would 
already have passed. From Fig. 6 one sees that at any 
given time AGASA would record only part of the total 
energy spectrum, since particles with = 10 EeV arrive 
at time r ^ lO^yrs, while particles with E = 100 EeV 
arrive at times t ~ 10® yrs. 

This implies that if Cen-A is the source of all UHE- 
CRs, it has to be a continuously emitting source, or, 
what amounts to the same, an intermittent source that 
emits on timescales tem, with quiescent periods of dura- 



tion At ^ min[Ar(i?)]. Here At{E) denotes the spread 
of time delays at energy E, and the condition ensures 
that the various contributions of the various bursts of 
particles largely overlap so as to produce a featureless 
energy spectrum at all times. Furthermore, since no 
low energy cut-off in the energy spectrum has been seen 
down to energies ~ 5 EeV, Cen-A must have been active 
for a time corresponding to the largest time delay pos- 
sible ~ r(5EeV) + AT(5EeV). For S ~ 1 /iC, we find 
that Cen-A must have been producing UHECRs inter- 
mittently for the past ~ 10 Gyrs. Whether this is realistic 
or not can hardly be constrained on theoretical grounds. 
The above constraints on the duration of the periods of 
activity and quiescence tcm and At read tem ^ 10^ yrs, 
and similarly for At, which are reasonable orders of mag- 
nitude for the evolution timescale of Cen-A [26] . The age 
of Centaurus A can be estimated from the time necessaxy 
for its jets to extend to their present size ~ 250 kpc, i.e. 
between ~ 10® yrs and a few Gyrs, depending on the de- 
celeration during extension, so that for our purposes this 
age is essentially unknown. Note that Centaurus A is a 
radio-galaxy with sub-relativistic jets, and without hot 
spots; the lobes are not particularly active, with a total 
bolometric luminosity ~ lO'^^ergs/s [26]. 

The constraints of Ref. [19] on the energy requirement 
must thus be reconsidered. In the stationary regime, 
the flux at Earth reads E^j{E) = E'^q{E)/{A'KfD{E)d, 
where q{E) is the injection spectrum at the source. One 
easily calculates that, assuming q{E) oc E~'^ , in order 
to produce the energy weighted flux measured by cosmic 
ray experiments ~ lO^**'^ eV^m~^s~"'^sr~^, one requires a 
UHECR emission power Puhecr — lO^^SZgdcenA ergs/s 
at the source. Note that this represents the average 
power; the actual power during the phase of activity is 
higher and reads 'PuHECR(iom + AT) /tem- The above 
"PuHECR is thus a strict lower limit, which is neverthe- 
less more optimistic than the ~ 10"^^ ergs/s obtained in 
Ref [19]. 

As an aside, we note that the above constraint on the 
age of Centaurus A can be generalized to any such source 
of UHECRs, namely that the time delay at the lowest en- 
ergies be smaller than the age of the source, and thus also 
the age of the Universe. The propagation time at 5 EeV 
reads, using Eq. (3): r ~ 2.8 Gyrs B-e^cenA- Imposing 
T < 14 Gyrs gives a general constraint between the dis- 
tance d to the source and the strength of the intervening 
magnetic fleld: 

Even if there are several sources contributing to the 
cosmic ray flux, this limit should hold imless the sources 
conspire to add their individual piecewise contribution 
in such a way as to form a featureless energy spectrum. 
When there are many sources the above constraint disap- 
pears, as the central limit theorem would guarantee that 
a featureless spectrum would be produced; this is no- 
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tably one of the peculiarities of the 7— ray burst model 
of UHECR origin [27]. 

IV. CONCLUSIONS 

Our detailed numerical simulations show that the 
model considered in Ref. [19], in which Centaurus A 
is the source of all observed UHECRs, is inconsistent 
with the data, at least for the magnetic field strength 
B ~ 0.3 /LtG put forward by these authors. We find that 
for a magnetic field strength B ~ 1 /iG, the predicted 
energy spectrum is in marginal agreement with that ob- 
served by AGASA. However the large deflection angle of 
the highest energy event (the Fly's Eye event) with re- 
spect to the line of sight to Ccn-A must be explained as 
a ~ 2(7 fluctuation. We also argued that this magnetic 
field intensity saturates the observational upper bounds 
from Faraday rotations and on X-ray emission from the 
ambicn gas. This model can be tested by improving 
these upper limits with current experiments. We also 
showed that in order to explain all UHECRs down to 
E ~ 5 EcV, Cen - A must have been producing UHECRs 
for the past ~ lOGyrs. All these facts are rather strong 
requirements on the source and on the intervening mag- 
netic flelds. We thus find anew a conclusion obtained in 
previous studies [10,15 17], namely it is much more likely 
that a few sources within ~ 10 Mpc from the Earth would 
produce the observed ultra-high energy cosmic ray flux, 
and that the ambient magnetic field strength in the Lo- 
cal Supcrcluster B ~ 0.1 — 0.5 /iC. Work is in progress to 
quantify the number of sources needed and the distance 
scale for various values of the magnetic field strength. 
Ongoing and future ultra-high energy cosmic ray experi- 
ments [5,6,21], by increasing the statistics at the highest 
energies, will soon provide much tighter bounds on the 
number of UHECR sources. 
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